Dynamically-Induced Frustration as a Route to a Quantum Spin Ice State in 
Tb 2 Ti 2 7 via Virtual Crystal Field Excitations and Quantum Many-Body Effects 
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The Tb2Ti207 pyrochlore magnetic material is attracting much attention for its spin liquid state, 
failing to develop long range order down to 50 mK despite a Curie- Weiss temperature #cw ~ 
— 14 K. In this paper we reinvestigate the theoretical description of this material by considering a 
quantum model of independent tetrahedra to describe its low temperature properties. The naturally- 
tuned proximity of this system near a Neel to spin ice phase boundary allows for a resurgence of 
quantum fluctuation effects that lead to an important renormalization of its effective low energy 
spin Hamiltonian. As a result, Tb2Ti207 is argued to be a quantum spin ice. We put forward an 
experimental test of this proposal using neutron scattering on a single crystal. 



Magnetic frustration arises when the lattice geometry 
prevents a system from finding its classical ground state 
energy by minimizing the energy between pairs of inter- 
acting magnetic moments (spins), pair by pair. Partic- 
ularly interesting are models of geometrically frustrated 
magnets where there exists a macroscopic number of clas- 
sical ground states not related by any global symme- 
try^. A prominent class of such systems are the spin 
ices where Ising spins reside on a three-dimensional py- 
rochlore lattice of corner-sharing tetrahedra 2 ^^. Because 
of their macroscopic number of quasi-degenerate low- 
energy states, spin ice materials possess an extensive low- 
temperature magnetic entropy^i similar to that found 
in the proton disorded phase of common water ice^. 

A current and exciting direction of research in frus- 
trated magnetism is the study of low energy effective 
Hamiltonians and gauge theorie o 9 ' 10 ' 11 to describe highly 
frustrated systems which, when ignoring quantum ef- 
fects, display an extensive classical ground state degen- 
eracy similarly to spin ices. Despite the seemingly broad 
conceptual context of gauge theory approaches, there 
have so far been few real frustrated quantum magnetic 
materials identified as potential candidates for the ex- 
otic behaviors proposed by these theories 1 ^. In this 
paper we argue that the paradoxical Tb 2 Ti 2 0y (TTO) 
pyrochlor o 12 i 13 ' 14 i 15 ' 16 i 17 ' 18 belongs to such a class of ma- 
terials. Specifically, we use a simple model to illustrate 
that the starting point of the above theories, namely a 
frustrated Ising spin ice Hamiltonian plus weak trans- 
verse terms, indeed constitutes the predominant part of 
the low energy effective Hamiltonian, H c s, of TTO. How- 
ever, as we show below, the microscopic mechanism that 
leads to the "dynamically-induced frustration" and the 
proposed spin ice assignment for TTO has heretofore es- 
caped scrutiny. We find that frustration and the spin- 
ice-like structure of i? c g dynamically emerge from vir- 
tual transitions to excited single-ion crystal field (CF) 
states and, most importantly, from quantum many-body 
effects. These transitions drastically modify the sym- 
metries of the many-body wave functions in the low en- 



ergy sector, leading to a significant renormalization of 
the longitudinal (Ising) part of H c g. This renormaliza- 
tion plays a crucial role for materials, such as TTO, that 
are naturally-tuned near the boundary between a Neel or- 
dered phase and the spin ice states. In particular, these 
transitions reposition TTO in the spin ice region of cou- 
pling parameter space. We are led to suggest that TTO 
is a novel quantum variant of the classical Ising spin ice 
materials studied so far— i2»2isL2j!. 

The main reason for the interest devoted to TTO lies 
in its failure to develop long-range order down to at least 
50 mK despite an antifcrromagnetic Curie- Weiss tem- 
perature, #cw ~ -14 Ki 3 -. Similarly to the Dy 2 Ti 2 07 
(DTO) and Ho 2 Ti 2 7 (HTO) spin ices, magnetic Tb 3+ 
in TTO possesses a single-ion CF Ising ground state 
doublet with wavefunctions "Fq ) where (^IJ 2 !^) are 
the only non-vanishing matrix elements of the J an- 
gular momentum operato r 19 i 20 . Monte Carlo simula- 
tions of a model with such classical Ising spins^ that 
can only point "in" or "out" of an elementary tetrahe- 
dron^ and interact via nearest-neighbor (nn) antiferro- 
magnetic exchang e) 21 ' 22 and long-range dipolar couplings 
predict, in dramatic contrast with the experimental find- 
ing o 12 i 14 ' 15 i 16 , a transition to a four sublattice Neel order 
at T c ~ 1.2 K§. A key difference between TTO and spin 
ices has so far not been carefully investigated: in spin 
ices, the excited CF states lie at an energ y 13 ' 19 several 
hundred times larger than the exchange and dipolar in- 
teractions and there is therefore little admixing between 
the CF states induced by the spin interactions. This is 
not necessarily the case for TTO where the first excited 
doublet lies at only A ~ 18.7 K above the ground Ising 
doublet 1 ^. It is therefore necessary to investigate how the 
-ffeff of TTO is affected by virtual quantum mechanical 
CF excitations. 

The Hamiltonian of TTO is taker>2i as H = H c{ + 
H + Hd- H c { is the single-ion CF Hamiltonian 1 ^, H e 



= JE 



is the nn exchange interaction and 
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where Ri is the position of atom i with total angular 
momentum J^. J7 is the nn exchange coupling with 
the convention here that J > is antiferromagnetic. 
V = (^ /47t)(5Mb) 2 /-Rnn i s the dipolar coupling, and 
g = 3/2 is the Lande factor for Tb 3+ . R nn = 3.59 A 
is the nn distance, giving T> « 0.0315 K — . To in- 
troduce the single-ion wavefunctions which become ad- 
mixed by the spin interactions, flint = fl + fld, we 
focus on the essential part of fl c f: its doublet ground 
states, "Jq ), and its lowest excited doublet states, 
at an energy A = 18.7 K above The excited 

states above A do not lead to qualitatively different new 
physics. Tb 3+ has orbital angular momentum L=3, spin 
S=3, and total angular momentum J = L + S with J=6. 
We express \^^) and \^f^) in terms of the eigenstates 
J = 6,toj) of J z within the fixed J = 6 manifold. Ex- 
act diagonalization of fl c f using the CF parameters taken 
from Ref. [H for Ho 2 Ti 2 7 but rescaled for Tb 2 Ti 2 7 
gives: |*±)= a 4 | ± 4} ± a 5 | T- 5) + a 2 | T 2} ± a t \ ± 1} 
and |*±> = ±/? B | ± 5) + /3 4 | ± 4> + /3 2 | ± 2) ± /? x | =F 1). 
With {a 4 ,a 5 ,a 2 ,ai}= (0.9402,-0.2958,0.1389,0.0960) 
and {A5,At,/3 2 ,/?i}= (0.9271,0.3117,-0.1917,0.0809). 

Exact diagonalization — single tetrahedron. Since the 
spin correlations in TTO never exceed a length scale 
much beyond a single tetrahedron^ 5 -, we first consider 
a simple model to describe TTO which consist of non- 
interacting tetrahedra. Such an approximation of inde- 
pendent tetrahedra explains semi-quantitatively several 
bulk properties of the classical Heisenberg pyrochlore an- 
tiferromagnet model 2 - . The same approximation also ac- 
counts qualitatively well for the spin-spin correlations^ 
which, by incorporating transverse spin fluctuation !] 22 ' 24 , 
captures the rough features of the neutron scattering of 
TT0 15 ». Our aim in using this approximation is to expose 
the general effects of virtual CF excitations on fl e ff . We 
ignore the long-range dependence of the dipole-dipole in- 
teractions in fldip since it is now well understood that it 
is the nn contribution of the dipolar interactions that pre- 
dominantly controls the transition from the Neel phase to 
the spin ice state^ 5 -^. Henceforth, we set V = 0.0315 K 
and fix the values of {a m } and {/3 m } to those listed 
above. We then treat J and A as independent tunable 
parameters in order to expose the essential physics at 
play near the Neel — spin ice boundary. 

Diagonalizing fl; nt for a single tetrahedron within the 
space of the 4 4 =256 CF states, we obtain the zero tem- 
perature J — A phase diagram shown in Fig. 1. For 
the classical Ising limit (1/A = 0), we recover the tran- 
sition between "Neel order" (all- in/all-out, two-fold de- 
generate, J > 52?) and a spin ice manifold (two-in/two- 
out, six fold degenerate, J < 52?) at J nn = Am, where 
J„n = i7K* + |J1* +)| 2 and D nD ee |2?|(*+|J-|vI/+)| 2 
Refs. '6.27}. A classical Ising model&i 3 - places TTO above 
the classical J nn = -Dnn boundary&22 (horizontal dashed 
line), i.e. in the Neel stated. However, for 1/A > 0, 
quantum fluctuations, due to the admixing of with 
) via flint, become increasingly important, as shown 
by the renormalized J7c(l/A) boundary in Fig. 1 (filled 
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FIG. 1: (Color online). J — A phase diagram of a single 
tetrahedron. TTO has J = 0.167 K and A = 18.7 K (cross 
symbol). Main panel: The boundary J c (l/A) (filled circles) 
separates a ground state singlet (J < J c ) from a ground state 
doublet {J > J c ). The open circles show the same bound- 
ary, but as predicted by exact diagonalization of H B s . Inset: 
neglecting transverse terms in H e s (A = 0), the filled squares 
show the renormalized classical sextet-doublet boundary set 
by the condition J"(j7, 1/A) = 0. The open squares show the 
incorrect sextet-doublet boundary predicted when Jff for pair 
in H e ff is obtained by ignoring contributions in PHRHP 
coming from (intermediate) excited states ) that belong 
to the two other Tb 3+ ions (k and I) on the tetrahedron. 



circles). This boundary separates quantum variants of 
the classical phases and TTO, with J = 0.167 K and 
A = 18.7 K (cross symbol in Fig. 1), is now deeply repo- 
sitioned in the singlet regime, i.e. is a quantum spin ice. 
The quantum spin ice state for J < J c is a singlet pre- 
dominently built of the fully symmetrized 6 two-in/two- 
out otherwise degenerate classical spin ice states whose 
degeneracy is lifted by quantum effects. The ground state 
also contains a small (of order 1/A) spectral weight con- 
tribution from the excited CF states. The ground singlet 
is accompanied by a low energy spectrum of 15 excited 
states that consists of three triplets and three doublets 
spanning an energy band SW w 0.5 K above the ground 
state, and which is separated by a gap of 16 K from 
other high energy states. Above the boundary J7c(l/A), 
the ground state is an all- in/all-out doublet similarly to 
the classical 1/A -> limit 3 .^. 

Single-tetrahedron model — effective Hamiltonian. To 
shed light on the virtual CF excitation channels lead- 
ing to the J" C (A) above, we construct an effective 
S = 2 anisotropic Hamiltonian, fl c ff. Using sec- 
ond order perturbation theory^, in 1/A, we have 
fleff = PHP + PHRHP, with P = Eal $ 0,«)<*0,a| 
and R = E/sl^Hflo-fl^r 1 ^!, where E = 
($0 )0! |flcf|$o,a) and Ep = ($ e ,/3|flcf |$ e ,/3>- Here {|$0,a)} 
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are the 2 



16 states constructed 



as direct prod- 



ucts of the non-interacting single ion |*jf ) CF doublet 
ground states of H c f. The |$ e ,/3) are the remaining 4 4 - 
16=240 states. We recast H e g in the form of an effective 
anisotropic S = i spin Hamiltonian in the individual 



local [111] spin a? basisSS: H cS = E { i, jhfi ,u J%°?°j, 
where [i, v are spin component indices, erf are Pauli ma- 
trices, and J^j are the effective anisotropic coupling con- 
stants. A constant energy term has been dropped from 
H c g, while the one-site J^o? terms get eliminated by the 
symmetry of a tetrahedron. Figure 1 shows the singlet- 
doublet boundary predicted by H c s (open circles). 

In order to expose the most important aspects 



of 



expose 
we write it 



as 
S 



E 
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with J;, 



Jij {J , 1/ A) and with the perturbation parameter A ul- 
timately set to A = 1. H c g contains transverse (non- 
Ising) - <^<^)of <rj terms where, for A = 18.7 K, 

the largest of the transverse J?? is approximately 50% 
of the Ising J zz coupling. To generate these terms via 
the nonvanishing matrix elements of Jf and J, between 
|$o,q) and I'&o./s) in PHRHP, it is important to re- 
tain more than the predominant ±4} and f3§\ ±5} 
components in \*5f^) and The term PHP cor- 

responds to the classical [111] Ising model with nn ex- 
change and dipolar interaction s^ 13 . It is by accident 
that J jD has a specific value such that PHP almost 
vanishes for TTO&2£, hence allowing an opportunity for 
the resurgence of quantum effects via PHRHP in if e ff- 
The contribution of PHRHP to J zz a z a z is ferromag- 
netic, and hence competes with the antiferromagnetic 
classical PHP Ising term and brings back frustration in 
TTO. Neglecting momentarily the quantum transverse 
terms (A — > 0), the change of sign of Jff(J , 1/A) con- 
trols the transition from a (spin ice) two- in/two-out sex- 
tet to an all- in/all-out doublet (filled squares in inset 
of Fig. 1). It is a key point of this paper that it is 
the renormalization of the Ising sector of the theory, 
Jif(J,fy -> JifiJ, 1 /^), caused by virtual CF exci- 
tations, that largely determines the ,7c (1/A) boundary 
(filled circles, main panel) and its upward movement with 
respect to the classical J Z 3 Z {J,Q) = boundary (hori- 
zontal dashed line). Specifically, compare the curve with 
filled squares in inset of Fig. 1 with the curve with filled 
circles in the main panel, and note the semi-quantitative 
agreement. 

It is important to note that the virtual excitation of 
an intervening (third) ion k, with angular momentum 
and H = H t k + Hkj in PHRHP, plays a crucial role in 
the renormalization of the classical Ising sector Jffafa z 
for pair (i, j) in H c s . Only by including this "third body" 
contribution do we get the correct trend for the "classi- 
cally renormalized" J ZZ (J,1/A) — boundary. Failure 
to do so gives an incorrect boundary decreasing with in- 
creasing 1 /A (curve with open squares in inset of Fig. 1). 
Hence, it is the quantum many body aspect of the full 
microscopic quantum H — H c f + H c + H^ of TTO that 
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FIG. 2: (Color online) Theoretical diffuse neutron scattering 
intensity I(q)/\F(q)\ 2 for a single tetrahedron with four CF 
states per Tb 3+ in (hhl) plane at 9 K. 



produces the interesting physics here, namely the correct 
renormalization of the Ising part of its iJ e ff from an un- 
frustrated system when 1/A = to that of a frustrated 
ferromagnetic nn spin ice model 3 -. The aforementioned 
transverse (quantum) part of H c s (A ^ 0) lifts the de- 
generacy of the (classical spin ice) sextet, giving a singlet 
ground state for an independent tetrahedron and 15 ex- 
cited states within SW « 0.5 K above the ground state. 

Diffuse neutron scattering. The result that iJ e fF for 
TTO is a nn ferromagnetic Ising spin ice model plus 
transverse terms might come as a surprise and be per- 
ceived as incompatible with neutron scattering measure- 
ment o 14 ' 15 ' 16 . Indeed, the neutron scattering pattern of 
TTQ 14 ' 15 i 16 is qualitatively very different from that of 
the HTO 22 and DTO ^ spin ices. In TTO, there 
is an intensity maximum at 002 in the (hhl) scatter- 
ing plane and a second broad maximum at 220. In 
spin ices, there are broad maxima at 003 and ||l 
— i 29 ' 30 . Therefore, the question is whether the above 
single tetrahedron model characterized by a ferromag- 
netic Ising sector in its H c g gives a diffuse neutron 
scattering pattern compatible with experimen t 14 i 15 i 16 . 
To address this question we compute the diffuse neu- 
tron scattering intensity, I(q), using standard formu- 

21 c( a .«; b ./3) 



lac 31 : I(q) ex |F(q)| 2 £. Q JS a p - ^|q|- 2 ]5, 



J difT 



where a, b are the sites on the tetrahedron, a, (3 are 
spin components and F(q) is the Tb 3+ form factor. 



Q(a,a;b,P) 
^diff 



-E n /k B T 



where the states {n, n'} are those whose energy E n 
falls within the experimental energy/frequency resolu- 
tion window of ~ 4.3 K over which the neutron scat- 
tering intensity is energy integrated™. These are, in- 
cidently, the same low energy states that span an en- 
ergy SW w 0.5 K above the ground state. Numeri- 
cal results for I(q)/\F(q)\ 2 at T = 9 K are shown in 
Fig. 2. One finds a good match in the symmetry of 
the theoretical pattern in Fig. 2 with the experimental 
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FIG. 3: (Color online) Theoretical diffuse neutron scattering 
intensity I(q)/\F(q)\ 2 along q = (hh2)] at 400 mK (left panel) 
and 40 mK (right panel) as a function of the nearest-neighbor 
antiferromagnetic exchange J for a single tetrahedron with 
four CF states per Tb 3+ . 

one in Fig. 6 of Ref. [H|. These results show that, de- 
spite a predominant ferromagnetic nn Ising component, 
H e g possesses sufficient low-energy transverse response 
(fluctuations) to account for the symmetry of the diffuse 
neutron scattering. These give, in particular, the inten- 
sity maximum at 002 that arise from spin fluctuations 
transverse to the local [111] Ising direction a 15 i 22 ' 24 . We 
propose that a scan of I(q)/\F(q)\ 2 along the (hh2) di- 
rection may be used to ascertain whether TTO is indeed 
in a quantum spin ice state at low temperatures. Figure 
3 shows that I(q)/\F(q)\ 2 along (hh2) has a broad maxi- 
mum at h=0 in the singlet/spin ice regime, J < 0.187 K 
for 1/A = 0.053 K _1 (see main panel, Fig. 1), while 
it has maxima at h — ±8(J) in the doublet regime for 
J > 0.187 K. The split Iih2 intensity line scan as a char- 
acterization of the underlying (antiferromagnet vs spin 
ice) ground state is sharper the lower the temperature. 

Going beyond the single-tetrahedron approximation, 
competing anisotropic interactions further than nearest 
neighbors are generated by virtual CF excitations. These 
lead to a 7J e ff which, at the classical level displays vari- 
ous long-range ordered spin ice states depending on J/V. 
Preliminary calculations find a long-range ordered Q = 



ferromagnetic (spin ice) state, similar to the one recently 
reported for Tb2Sn2 07 32 , which competes with the previ- 
ously identified Q = Neel order— and the Q = 001 long 
range ordered spin ice of the dipolar spin ice mode l 25 ' 27 . 

In conclusion, we have used a simple model of non- 
interacting tetrahedra to describe the low temperature 
properties of of the Tb2Ti207 magnetic pyrochlore ma- 
terial. The present work identifies a new mechanism 
for dynamically-induced frustration in a physical sys- 
tem which proceeds via crystal field (CF) excitations 
and quantum many body effects. More specifically, we 
uncovered that interaction-induced fluctuations among 
otherwise non-interacting single-ion CF states lead to 
a renormalization of the low-energy effective theory of 
Tb2Ti207 from that of an unfrustrated [111] pyrochlore 
Ising antiferromagnet 3 ' 6 to a frustrated nearest-neighbor 
spin ice model 3 ' 4 . The remaining transverse fluctuations 
lift the classical ice-like degeneracy and, at the single- 
tetrahedron level, the system is in a quantum mechan- 
ically fluctuating spin ice state, or resonating spin ice. 
The effects discussed here are likely responsible for some 
of the subtleties underlying the failure of this material to 
order at a temperature scale of 1 K— i 13 i 18 i 22 . Whether the 
true quantum mechanical ground state of the full lattice 
model of Tb2Ti207 is a semi-classical long-range ordered 
state with finite quantum spin fluctuations 33 , or a more 
exotic quantum ground state d 10 ' 11 , is a challenging but 
very exciting problem for future studies. 
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